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UNCEASSIFIED 


Design  and  development  of  laser  Doppler  velocimetry  based  on 
DSP  technique  for  blood  flow  measurement 

Yuan-Hsiang  Lin,  Yan-Yu  Chen,  Yuan  Ho,  Yi-Lin  Li,  and  Gwo-Jen  Jan* 

Department  of  Electrical  Engineering,  National  Taiwan  University 

ABSTRACT 

A graphic  user  interface  and  real-time  laser  Doppler  velocimeter  (LDV)  based  on  tlie  digital  signal  processor  (DSP)  had 
been  designed  and  developed.  The  hardware  setup  included  tlie  Michelson  interferrometer  optics,  photo-detector,  current  to 
voltage  converter,  AC  amplifier  and  filtering  circuits,  as  well  as  a DSP  module.  Tlie  software  system  on  DSP  module  was 
also  developed  to  access  data  and  to  perfonn  the  moment  weighting  algoritluns.  In  addition,  tlie  processed  data  was 
transmitted  to  the  personal  computer  and  advanced  analysis  could  be  achieved.  Tlie  velocity  measurement  using  developed 
LDV  device  was  calibrated  by  a mirror  moimted  on  a linear  vibrator.  The  outcomes  presented  high  linearity  and  good 
accuracy.  In  vitro  experiment  employing  tliis  LDV  system  was  also  carried  out.  The  results  showed  that  the  developed  LDV 
instniment  offered  a flexible  tool  to  investigate  the  blood  flow  of  microcirculation  system. 

Keywords:  Laser  Doppler,  Blood  flow.  Digital  Signal  Processor  (DSP) 

1.  INTRODUCTION 

In  recent  years,  scientists  have  devoted  to  investigating  how  tlie  circulation  diseases  influence  microvascular  flow'.  Some 
relations  were  found  and  early-stage  detection  of  these  diseases  can  be  acliieved  by  continuous  monitoring  of 
microcirculation.  Among  various  measuring  techniques,  laser  Doppler  velocimetry  (LDV)  has  been  adopted  most  widely. 
Since  first  introduced  by  M.D.  Stem^  in  1975,  LDV  is  now  well  established  and  extensively  used  in  skin  perfusion 
measmements.  Other  applications  to  monitor  choroidal  flow^,  renal  arterioles''  and  microcirculation  in  gums  were  also 
developed.  Because  tlie  tecluiique  is  uon-invasive  and  liighly  spatial  sensitive,  so  far  it  can  liardly  be  replaced  by  any  other 
methods  in  physiology  measurement. 

However,  though  LDV  provides  such  advantages,  there  are  still  some  problems  to  be  solved^.  Since  tlie  scattering 
mechanism  is  too  complicated  to  be  determined  in  tissue®,  various  algoritluns’  * are  developed  to  obtain  indicators  for  fluid 
velocity.  Applying  different  algorithms,  the  results  may  somewhat  differ.  In  addition,  tlie  processing  bandwidtli®  is  also 
involved.  Adopting  low  cut-off  frequency  may  lead  to  under-estimation  of  velocity  while  using  liigh  cut-off  frequency  limits 
tlie  frequency  resolution.  Besides,  tlie  computed  result  is  closely  related  to  tlie  penetration  deptli'",  which  is  decided  by  the 
laser  used".  With  different  penetration  depths,  different  volumes  are  sampled;  and  it’s  hard  to  say  which  interprets  the  fluid 
velocity  better.  More  tlian  tliat,  absolute  in  vivo  calibrations  of  LDV  systems  are  difficult  to  acliieve  due  to  tlie  large 
variation  in  cutaneous  stnicture”. 

Therefore,  in  this  study  we  developed  a digital  signal  processor  (DSP)  based  LDV  system  and  a simplified  physical  model 
was  also  implemented  for  the  in  vitro  experiment.  By  applying  the  DSP  framework,  we  evaluated  the  linearity  of  different 
algorithms,  and  determined  a linear  indicator  for  the  fluid  velocity. 
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2.  THEORY  AND  COMPUTATION  ALGORITHMS 

Wlien  a photon  is  scattered  by  a moving  particle,  the  light  frequency  is  slightly  shifted.  The  Doppler  frequency-shift  Af  is 
given  by’^ 

Af  = (|k||v|/ S' )sin'/2  0 cos  a (1) 

where  k is  tlie  incoming  wave  vector  ( k = 2 ;r  / X , A is  tlie  wavelength  ),  v is  tlie  velocity  of  the  moving  particle,  6 is  tlie 
scattering  angle,  and  a is  the  angle  between  Ak  and  v. 

Applying  Doppler  eSect  on  skin  perfusion  measurement,  tlie  situation  becomes  much  more  difficult  due  to  the  complicated 
scattering  mechanism  in  tissue.  Tlie  intensity  modulation  of  the  detected  signal  shows  an  approximately  exponential  decay 
in  tlie  spectra.  Boimer  and  Nossal’s''’  results  showed  tliat  tlie  first  ( weighted  ) moment  (&i)  of  the  spectral  power  density 
S(6j)  of  the  detector  signal  is  linearly  proportional  to  die  average  velocity  (v)  of  all  moving  particles  in  the  optical  sampling 
volume.  The  first  weighted  moment  (FWM)  is  expressed  as 

(v)  ~ (<y)  = Ml  / Mo,  with  Mn  = ^ co”S(w)da)-  (2) 

The  numerator  Ml  is  the  first  moment  of  the  spectral  power  density  S(a)).  It  is  assumed  to  be  proportional  to  the  blood  flow 
rate  in  the  optical  sampling  volume.  The  denominator  Mo  is  proportional  to  the  total  amount  of  Doppler-shifted  light. 

The  second  weighted  moment  (SWM)  is  expressed  as 

<v> ~ M /«,.  (3) 

J*  S(co)d(o 

Square  root  of  <v^)  represents  die  RMS  value  of  blood  velocity,  and  can  be  regarded  as  anodier  representative  of  the  blood 
velocity. 


After  the  calculation  of  the  first  and  second  weighted  moments,  the  results  can  be  derived,  in  die  unit  of  frequency  (Hz) 
instead  of  velocity  (m/s).  The  calculated  value  is  proportional  to  the  blood  velocity,  and  a conversion  factor  should  be  added 
to  compute  die  fluid  velocity.  However,  die  conversion  factor  varies  between  individuals  and  even  between  different  sites 
on  an  individual.  Thus  we  adopted  an  linear  indicator  to  represent  the  fluid  velocity  instead  of  the  absolute  value. 

3.  INSTRUMENTATION 


3.1.  Hardware  Setup 

Tlie  block  diagram  of  the  LDV  hardware  system  is  shown  in  Fig.  1.  It  includes  an  optical  system,  analog  signal  conditioning 
circuits,  a DSP  module  and  a graphic  user  interface  (GUI).  Tlie  optical  system  contains  a 5mW  laser  diode  (650nm)  and/or 
He-Ne  laser  (632. 8mn),  a test  object,  and  an  optical  detector.  It  detects  the  optical  heterodyne  and  converts  tlie  signal  into 
electrical  domain.  In  analog  module,  tlie  preamplifier  initially  amplifies  tlie  weak  signal  and  transfonns  photocurrent  to 
voltage  (gain  is  200  V/mA).  The  AC  amplifier  further  amplifies  tlie  AC  signal  and  the  gain  is  depending  on  the  light 
intensity  (max.  gain  is  200).  The  total  bandwidth  of  the  analog  module  is  50kHz. 

Tlie  DSP  module  consists  of  a 16  bit,  200  kHz  sampling  rate  ADC  for  A/D  conversion,  a 32k-byte  EEPROM  for  program 
storage,  liigh  speed  SRAMs  up  to  128k-word  for  DSP  data  buffer,  and  a TMS320C31  DSP  as  well  as  an  ADS7843 
touch-screen  controller.  In  tliis  module,  a 2048-points  Fast  Fourier  Transform  (FFT)  is  performed  on  DSP.  Tlie  sampling 
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rate  is  200kHz  and  tlie  frequency  resolution  is  97.7Hz.  According  to  equation  (1),  in  simulation  experiment  (A,=632.8nm) 
tlie  maximal  measurable  velocity  is  15.82imn/s  and  velocity  resolution  is  30.9pm/s  in  tliis  system  tliat  are  depended  on  6 
and  a . The  GUI  consists  of  a 240x  128  dots  graphic  LCD  and  a touch  screen  ptuiel.  Tlie  touch  screen  allows  a user  to  enter 
a command  simply  by  touching  a location  on  tlie  panel.  The  two  devices  together  make  the  user  interface  more  friendly  and 
easy  to  use.  In  tlie  other  hand,  tlie  DSP  can  transfer  tlie  results  to  PC  in  real-time  tlirough  RS-232  interface,  and  data 
analysis  and  storage  can  be  further  acliieved.  It  lias  flexibility  to  modify  the  computing  algoritlim  of  velocity  of  DSP 
through  the  control  panel  of  PC. 
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Fig.  1 The  block  diagram  of  the  LDV  hardware  system 


3.2.  Software  Development 


The  main  flowchart  of  our  program  is  shown  in  Fig.  2.  After  pow'er  on,  program  starts  at  initializing  all  control  registers  of 
DSP  and  LCD.  Then  program  goes  into  tlie  major  loop,  wliich  makes  tlie  DSP  wait  for  sample  data  to  fill  an  input  buffer, 
perform  FFT  on  the  data  to  obtain  the  power  spectrum,  calculate  blood  velocity  based  on  tlie  spectrum,  then  display  it  and 
carry  out  any  user  input.  For  the  requirement  of  real-time  processing,  all  source  codes  were  written  in  C3x  assembly 
language  instead  of  other  high-level  languages. 


The  FFT  used  in  the  program  is  a 2048-point,  decimation  in  time  algorithm.  Powered  by  the  bit-reverse  addressing  mode 
and  parallel  instnictions  of  the  C3x  DSP,  the  FFT  can  be  calculated  in  a very  short  period  of  time.  With  a sampling  rate  of 
200kHz,  the  frequency  resolution  df  is 


look 

2048 


91.1Hz . 


(4) 
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Sampling  of  signal  goes  continuously,  and  every  5psec  DSP  will  receive  a new  sample  data  from  ADC.  To  prevent  new 
data  from  mixing  with  data  tliat  are  currently  under  processing,  a three  buffer  rotating  teclmique  is  used.  Input  buffer 
accumulates  new  data  coming  from  ADC.  Tlie  data  in  FTT  buffer  are  tliose  currently  under  processing  (FFT  or  calculation 
of  blood  velocity),  and  output  buffer  stores  previously  processed  data,  which  is  used  for  the  spectrum  display. 


Fig.  2 The  flowchart  of  the  main  program 

When  sample  data  fills  the  input  buffer,  a buffer  rotating  occurs.  The  input  buffer  becomes  FFT  buffer  so  tliat  its  contents 
are  ready  for  FFT  calculation.  FFT  buffer  becomes  output  buffer  and  tlie  spectrum  stored  in  it  will  be  displayed  immediately. 
Output  buffer  becomes  input  buffer  and  will  start  accepting  new  sample  data. 

Tlie  calculation  of  blood  velocity  involves  tliree  algoritluns  -FWM,  SWM,  and  the  peak  value  of  the  frequency  (PVF) 
metliods.  PVF  was  devised  by  us  and  is  simply  tlie  selection  of  the  frequency  that  lias  maximum  power  in  the  power 
spectrum.  Since  tlie  power  spectrum  stored  in  memory  is  actually  in  discrete  fonn,  tlie  integration  of  equation  (2)  must  be 
approximated  by  the  following  summation: 

«2 

^ndf.S[n\ 

FWM  = . (5) 

Z-S’W 

n-n\ 

Where  df  is  the  same  as  in  equation  (4)  and  S[n]  is  tlie  power  spectrum  after  FFT  calculation.  It  is  obvious  tliat 
frequencies  lower  than  n\  ■ df  and  tliose  liigher  tlian  n2  ■ df  will  be  excluded  from  the  calculation  of  blood  velocity. 
Thus  by  choosing  the  value  of  n\  and  n2 , a tliresholding  mechanism  is  easily  implemented,  witli  n\  ■ df  and 
n2  ■ df  being  the  lower  and  upper  tlireshold.  Proper  use  of  tliis  tliresholding  meclianism  can  eliminate  unwanted 
frequency  components  like  high  frequency  noises  or  the  large  DC  component  in  the  signal  spectrum. 

Similarly  the  calculation  of  SWM  uses  the  following  summation  form  instead  of  the  integration  form  in  equation  (3). 


201 


SWM  = Sqrt 


(6) 


f »2  \ 

^{ndfy-S[n] 

n=n\ 

n2 


V n=nl  J 


Tile  PVF  algoritlun  was  iinplemented  using  a simple  compare-and-update  algorithm  for  finding  the  frequency  that  has 
largest  power  in  the  signal  spectrum. 

4.  EXPERIMENTS 


4.1.  Simulation  Experiment 


FWM,  SWM,  and  PVF  metliods  are  designed  to  calculate  tlie  velocity  of  a moving  object.  The  simulation  experiment  is  to 
check  tlie  computation  of  tlie  algoritluns  and  to  calibrate  the  LDV  device.  Tlie  experimental  setup  is  shown  in  Fig.  3.  A 
632. 8mn  He-Ne  laser  provides  a light  beam  wliich  pass  tlirough  tlie  beam  splitter  (50:50)  and  splits  into  two.  One  beam  liits 
tlie  moving  mirror  and  is  Doppler  frequency  shifted.  The  other  hits  a stationary  mirror  and  keeps  its  frequency  non-shifted. 
The  reflected  beams  are  collected  by  an  optical  fiber  and  picked  up  by  a detector. 


Fig.  3 The  block  diagram  of  measurement  system  setup  for  simulation  experiment 

Tlie  iiioving  mirror  is  mounted  on  the  front  surface  of  a speaker  center.  A function  generator  outputs  a triangular  wave  to  tlie 
speaker  and  vibrates  the  mirror.  In  the  forward  and  backward  motion  of  tlie  speaker,  the  mirror  on  its  center  would  have  a 
constant  velocity.  The  relationsliip  between  applied  voltage  and  speaker  displacement  has  be  measured.  The  linear 
relationship  from  the  experiment  ineasureinent  is  shown  in  Fig.  4. 

For  the  forward  motion,  die  speaker  kept  pusliing  the  mirror  forward  in  velocity  V , which  can  be  calculated  from 
frequency  and  peak-to-peak  voltage  of  the  triangular  wave.  For  convenience  the  backward  motion  is  neglected.  The  trigger 
signal  from  the  fiuictioti  generator  is  connected  to  the  LDV,  which  could  be  set  to  trigger  mode  so  that  it  functioned  only  for 
the  forward  motion  of  the  speaker,  be  sure  the  measured  velocity  is  constant.  If  the  shifted  frequency  hasn’t  be  triggered,  tlie 
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peak  would  vibration  in  tlie  spectrum.  Tliis  fimction  also  could  be  used  to  investigate  tlie  sequence  relation  of  ECG  signal 
and  the  blood  flow  in  microcirculation. 

Applying  Doppler  frequency  shift  equation  (1)  we  get 

^if{kHz)=^3.\6\xv{mmls)  (7) 

where  A/  is  Doppler  frequency  shift,  v is  Mirror  velocity. 

y = 0.2712x  + 1.3181 
R’  = 0.9999 
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Fig.  4 Displacement  calibration  curve  of  a vibrating  mirror 

Note  tliat  tliis  conversion  factor  between  frequency  and  velocity  is  valid  only  in  tliis  simulation,  and  cannot  be  applied  to 
any  other  situations.  By  changing  the  peak-to-peak  voltage  of  the  triangular  wave,  tlie  velocity  of  the  mirror  changed  as  well. 
For  each  velocity  the  outputs  of  the  tliree  algoritlims  (FWM,  SWM,  and  PVF)  were  recorded.  Tlie  results  are  discussed  in 
tlie  following  section. 

4.2.In  Vitro  Experiment 

In  order  to  evaluate  our  system,  we  developed  a simplified  close-loop  model,  in  wliich  tlie  linearity  of  various  algorithms 
can  be  tested.  The  schematic  diagram  of  the  setup  is  shown  in  Fig.  5.  In  tliis  model,  micro-spheres  were  added  into  the  water 
to  serve  as  scattering  particles  (mass  concentration  = 0. 1 %),  and  a DC  motor  witli  a precise  controller  was  also  employed  to 
pump  the  solution.  The  photo-detector  of  the  LDV  picked  up  the  light  signals  scattered  either  by  tlie  moving  particles  or  tlie 
static  tube.  The  Doppler  shifted  and  non-sliifted  parts  of  the  scattered  beams  heterodyned  at  the  photodiode  and  produced 
corresponding  photocurrent.  After  being  amplified,  tlie  analog  signals  were  translated  into  digital  signals  and  a FFT 
processing  was  performed.  Deriving  tliis  iiifonnation  in  Uie  frequency  domain,  we  applied  different  algoritluns  to  search  for 
a linear  transformation  between  the  frequency  spectnun  and  tlie  fluid  velocity.  The  processed  data  were  transmitted  to  tlie 
PC,  and  advanced  analysis  could  be  taken. 
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Fig.  5 Experimental  setup  for  in  vitro  measurement 


5.  RESULTS  AND  DISCUSSION 

5.1.  Simulation  Experiment 

The  design  value  of  the  mirror  velocity  can  be  calculated  from  the  displacement  measurement  and  vibrating  period.  Tlie 
measured  values  using  the  FWM,  SWM  algorithms  and  PVF  method  are  originally  in  terms  of  frequency  unit  (Hz)  and  are 
converted  to  the  unit  of  velocity  (mm/s)  using  equation  (7).  The  comparison  between  designed  value  and  measured  value 
was  displayed  in  Fig.  6 (A),  (B),  and  (C).  As  the  figtue  shows,  all  tlmee  algoritluns  gave  values  lower  than  tlie  design  value. 
This  is  because  the  actual  mirror  velocity  is  lower  than  the  design  value.  Since  speaker  displacement  versus  applied  voltage 
is  calibrated  under  DC  condition,  when  an  AC  signal  is  applied  to  the  speaker,  tlie  speaker  will  exliibit  a frequency  response 
and  its  displacement  will  tend  to  be  smaller  than  in  DC.  Tlierefore  the  real  mirror  velocity  is  smaller  tlian  tlie  designed  value 
and  should  be  very  close  to  the  measiued  values  in  Fig.  6.  If  tlie  frequency  of  AC  bias  is  decreased,  the  difference  be 
reduced. 


01234S6  0123456  0123456 

Mirror  Velocity  (rmVs)  Mirror  Velocity  (rm/s)  Mirror  Velocity  (rmys) 


Fig.  6 Measured  velocity  versus  mirror  velocity  for  the  (A)  FWM  (B)  SWM  algorithm  and  (C)  PVF  method 
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However,  due  to  die  indetermination  of  die  conversion  factor  described  in  section  2,  die  important  parameter  is  not  the 
absolute  value  of  an  algoridim’s  output.  It  is  the  linearity  with  which  algorithm  responses  to  the  simulator’s  velocity,  tliat  is, 
how  straight  the  curves  in  Fig.  6 are.  Table  1 shows  the  coefficient  of  determination  of  the  four  sets  of  data.  is  an 
indication  of  linearity  of  a set  of  data,  with  R*  = 1 means  perfectly  linear  and  R^  = 0 means  totally  unrelated.  It  can  be 
observed  from  Table  1 that  all  three  algorithms  yield  very  good  linearity  with  the  mirror  velocity. 


Table  1 Coefficient  of  determination 


Algorithm 

Coefficient  of 
Determination 

First  Weighted  Moment 

0.9991 

Second  Weighted  Moment 

0.9988 

Peak  Value  Frequency 

0.9991 

5.2.  In  Vitro  Experiment 

The  bandwidth  is  cut-off  on  800Hz  and  line  width  is  IHz  in  our  experiment  setup.  Three  groups  of  spectral  response 
corresponding  to  different  rotation  rates  are  shown  in  Fig.  7.  Widi  an  increase  in  fluid  velocity,  die  fractions  of  the  Fourier 
components  are  reduced  at  low  frequencies  and  increase  at  high  frequencies. 


Frequency  (Hz) 

Fig.  7 Three  groups  of  spectra  corresponding  to  different  velocities  processed  by  moving  average  (10  periods) 

Fig.  8 shows  the  computed  velocities  versus  the  DC  motor’s  rotation  rate.  For  each  rotation  rate,  20  samples  were  taken  and 
averaged.  According  to  the  highly  correlated  relation  between  die  rotation  rates  and  the  computed  values  employing  FWM 
algoritlim  (R^=0.9922),  the  derived  value  might  be  a reliable  indicator  to  evaluate  the  relative  fluid  velocity.  An  indicator 
applying  SWM  is  also  provided.  Though  quite  linear  at  high  velocity  measurement,  it  seemed  to  deviate  at  low  velocity. 
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(A) 


(B) 


Fig.  8 (A)  Relationship  between  FWM  (in  arbitrary  unit),  (B)  SWM  (in  arbitrary  unit)  and  Rotation  Per  Minute  (RPM)  of  DC  motor 

6.  CONCLUSIONS 

We  have  developed  a portable,  real-time,  and  turn-key  LDV  system  based  on  DSP  teclmiqiies.  Two  most  frequently  used 
algoritluns  are  adopted  to  test  their  linearity.  Our  outcome  showed  that  both  algoritluns  provide  good  indicators  for  fluid 
velocity  measurement,  and  FWM  seems  to  be  a little  better  tlian  SWM.  Since  different  designs  of  simulation  model  and 
probe  may  also  influence  the  corresponding  results,  more  experiments  should  be  taken  to  further  evaluate  the  algorithms. 

Now  the  in  vivo  experiment  is  taken  on  a cow.  After  installing  a left  ventricular  assistant  device  on  the  cow,  we  use  the 
LDV  system  to  monitor  its  microvasculatine.  Information  collected  is  used  to  evaluate  tlie  target’s  recovery,  and  helps  tlie 
doctors  to  adjust  the  process  of  surgery. 
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